Dynamics of semifluxons in Nb long Josephson 0-n junctions 
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We propose, implement and test experimentally long Josephson 0-7r junctions fabricated using conventional 
Nb-AlOz-Nb technology. We show that using a pair of current injectors, one can create an arbitrary discontinuity 
of the Josephson phase and in particular a 7r-discontinuity, just like in d-wave/s-wave or in d-wave/d-wave 
junctions, and study fractional Josephson vortices which spontaneously appear. Moreover, using such junctions, 
we can investigate the dynamics of the fractional vortices — a domain which is not yet available for natural 
0-7r -junctions due to their inherently high damping. We observe half-integer zero-field steps which appear on 
the current-voltage characteristics due to hopping of semifluxons. 
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Theoretical and experimental investigation of Joseph- 
son junctions (JJs) made of unconventional superconductors 
showed that one can get so-called tt Josephson junctions! 1], 
for which the first Josephson relation has the form I s — 
—I c sm(f> = 7 c sin(0 + tt), instead of I s — I c s'm(j>. There 
can be different reasons for this. For example, in JJs formed 
using anisotropic superconductors with d-wave order param- 
eter symmetry the additional phase shift of tt happens when 
charge carries enter the negative lobe of the order parameter 
and exit from the positive lobe. Experimental realizations of 
7r -junctions include junctions formed at the boundary between 
two crystalline films of cuprate superconductors with different 
orientations ramp junctions between <i-wave and s-wave 
superconductors|301 or JJs with ferromagnetic barrier 010]. 

If one considers a ID long Josephson junction (LJJ) made 
of alternating and tt parts, half-integer flux quanta (so- 
called semifluxons|7, 8]) may spontaneously form at the joints 
between and tt parts. Semifluxons were observed using 
SQUID microscopy in different types of 0-tt- 

ljjsHhei. 

Semifluxons are very interesting objects which have not 
been studied in detail, especially experimentally. They are 
very different from integer fluxons. Fluxons are solitons, but 
semifluxons are not — they are spontaneously formed and 
pinned at the 0-7r-boundary. Consequently a semifluxon of ei- 
ther polarity represents the ground state of the system, while 
a fluxon moving in the LJJ always represents an excited state. 
This makes semifluxons attractive e.g. for information storage 
and processing in classical and quantum regimes. In spite of 
pinning, semifluxons are able to hop from one 0-tt -boundary 
to the other provided that they are not very far from each other. 
A chain of semifluxons usually exists in the antiferromagnetic 
ground state but can be "polarized" by means of an applied dc 
bias current fTHl . A semifluxon has an eigen-mode oscillation 
frequency, so a ID array of semifluxons has an energy band 
similar to the one in real crystals. By controlling the spac- 
ing between semifluxons (length of and tt -pieces) one can 
control the properties of such an artificial ID crystal. 

In this letter we propose and implement LJJs based on con- 



ventional superconductors which allow one to study arbitrary 
fractional vortices 1 25] and, in particular, their dynamics. Up 
to now this is not possible using natural 0-tt -junctions which 
typically have a Stewart-McCumber parameter j3 c ~ 5 or even 
smaller. Having much lower damping, we succeeded to ob- 
serve half-integer zero-field steps in the current-voltage char- 
acteristics which appears due to a non-trivial effect — hopping 
of semifluxons between neighboring 0-7r-boundaries. Low 
damping is also important in terms of future observations of 
macroscopic quantum effects involving fractional quanta| 12]. 

The dynamics of the Josephson phase in LJJs consisting of 
alternating and tt parts can be described by the ID perturbed 
sine-Gordon equation! 7] 

<pxx - 4>tt - sin (j> = a4> t -7 -9xx(%), (1) 

where c/>(x, t) is the Josephson phase and subscripts x and t 
denote the derivatives with respect to coordinate x and time t. 
In Eq. Q the spatial coordinate is normalized to the Joseph- 
son penetration depth A,/ and the time is normalized to the 
inverse plasma frequency w^ 1 ; a = 1/ \fjJ^ is the dimension- 
less damping; 7 = j/j c is the external bias current density 
normalized to the critical current density of the junction. The 
function 6(x) is a step function which is 7r-discontinuous at all 
points where and tt parts join. According to Eq. Q, <j){x) is 
also 7r-discontinuous at the same points as 9(x). We call these 
points phase discontinuity points. 

To describe 0-7r-LJJs one can use directly the equation with 
alternating critical current density ibi IsL 1 1 3l 1 1 411 

\i xx - ptt ± sin p = ap t - 7(1) (2) 

written for the continuous phase n{x,t) = 4>{x,t) — 0(x). 
Equations Q and (|2j are equivalent! 71 Ii26ll. 

Is it possible to create a 0-7T-LJJ (LJJ with 7r-discontinuity) 
using a conventional LJJ? The fact that a junction is 0-tt is 
contained in the 9 xx (x) term. If 6{x) is a set of steps, then 
xx {x) is a set of — S(x) / x singularities as shown in Fig.Qa)- 
(c). If we do not have initially the 9 XX term in Eq. ([0, we 
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FIG. 1: 8(x), 8 x (x) and 9 xx (x) and its approximation *yo(x) by two 
rectangular pulses of the width Aw and amplitude ty/Aw 2 . 



can substitute its effect by introducing a properly shaped bias 
current jg (x) = 6 XX (x) in addition to the main bias current 
7. To mimic one ^-discontinuity the 70 (x) profile should be 
the one shown in Fig.^c). For the sake of practical treatment, 
one can approximate it by two rectangular pulses, as shown in 
Fig. Qd). Such a bias profile can be created by two current 
injectors of width Aw situated as close as possible to each 
other. The current of the amplitude 7j n j = ttI c Xj/ Aw flowing 
from one injector to another should create a 7r-discontinuity of 
the Josephson phase. Thus one can create an artificial 0-7T-LJJ 
provided that the width of the whole injector construction is 
well below A,/. Note that, passing different currents, one can 
create arbitrary ^-discontinuities instead of a 7r-discontinuity, 
and study arbitrary fractional vortices, if they appear ll IB . This 
concept is a generalization of the idea to use a pair of injectors 
to insert a fluxon (27r-discontinuity) in an annular LJJ Ilal . 

In conventional LJJs the phase <fr = /1 is continuous while 
^inj = 0. When we increase 2i n j the phase (j> develops a jump 
k oc Jjnj at x — Xi n j. In practice, instead of a jump we have a 
rapid increase of the phase from <j) to (f> + k over a small but 
finite distance. Physically this means that by passing a rather 
large current through the piece of the top electrode between 
two injectors we "twist" the phase <j> by k over this small dis- 
tance, i.e., we insert a magnetic flux equal to k&q/(2it) in a 
small distance between injectors. This flux has a character- 
istic size of the order of the distance between injectors or of 
their width. The junction may react on the appearance of such 
a discontinuity by forming a solution 4>(x) [or fJ,(x)] which 
changes on the characteristic length Aj outside the injector 
area, and corresponds, e.g., to the formation of a fractional 
Josephson vortex, as discussed below. 

In the experiment the injectors are never ideal due to tech- 
nological constraints and one has to calibrate them to know 
the value of current needed to get, e.g., a ^-discontinuity. We 
do this by measuring the dependence of the critical current I c 



of the LJJ on the current through the injectors I- m j, and com- 
paring the obtained I c (Ji n j ) dependence with the I c (k) depen- 
dence calculated theoretically. As a result we get /j n j / k. 

The 7 C (k) curve (p/ c is the maximum normalized supercur- 
rent) numerically calculated for L = 2Xj (like in our sam- 
ples) is shown in Fig. [2] Details on the calculation of j c (k) 
for different L will be given elsewhere | la]. It is important to 
note that j c (k) has maxima at k = 2im and cusp-like minima 
(possibly with hysteresis) at n = 7r(2n+l). The value of -f c at 
the minimum depends on the junction length and varies from 




+2k 



phase discontinuity k 



FIG. 2: Normalized critical current 7 C vs. the Josephson phase dis- 
continuity k oc Jinj, numerically calculated for L = 2Aj. 



for a short junction to 2/ir for an infinitely long one| 16]. 

The samples were fabricated! 27] using standard Nb-AlOa,- 
Nb technology with low critical current density j c w 
100 A/cm 2 to have Aj ~ 30 fim as large as possible in com- 
parison with the injector size. Data reported here were ob- 
tained from the sample shown in Fig.[5]at T w 5 K. 

The measured I-V characteristics (IVCs) and I c (H) depen- 
dences (presented below) ensure good sample quality and the 
absence of trapped magnetic flux. 

To calibrate the injectors, we measured the dependence 
J c (J; n j) at H = shown in Fig. 0] Obviously, the depen- 
dence is similar to Fig. [5] and the minimum I c is about 10 % 
of the maximum I c , in agreement with numerical calculations 
for L/Xj = 2. From Fig. |4] we find that k = +ir corre- 



sponds to 



-13.5 mA and k 



-7r corresponds to 



w —12.8mA. The pattern is slightly shifted along the 
axis probably due to self-field effects. Thus, the maxi- 



mum I r is reached at Is, 



-0.3 mA. 



The value of I c at the maxima slowly decreases with in- 
creasing |/inj|- This effect can be reproduced in numerical 
simulations using the model Q for injectors with finite Ax 
and Aw like in our samples. 

One of the specific features of natural 0-7r -junctions is the 
minimum in the I C {H) dependence at H = o jill . It was 
calculated! 13] that when the 0-7T -junction is short (L <C Xj) 
the minimum is very deep and I c reaches zero at H = 0. As L 
increases, the minimum gets more shallow and at L ~ 10Aj 




FIG. 3: Optical microscope picture of the LJJ with two injectors. 
The width of the LJJ is w = 5 /im. 
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FIG. 4: Experimentally measured dependence of I c (Ii n j) (full cir- 
cles) and the maximum current of the semi-integer zero field steps 
^m(-Zinj) (open circles) at H = 0. The dependences are perfectly 
mirror symmetric relative to the Ji n j axis. 

almost disappears. In our artificial 0-7T-LJJ with L « 2Xj 
the minimum should be strongly developed. We have mea- 
sured and compared I C (H) at k = (no discontinuity) and at 
I C (H) at k = ±7r (I inj = +13.5 mA and I inj = -12.8 mA). 
The results are presented in Fig.[3]and are in good agreement 
with the prediction 1 13]. One can notice that at each value of 
Jinj corresponding to a 0-7T-LJJ, the left and right lobes of the 
I C (H) curve look somewhat asymmetric with the period of 
oscillations being different for H > and H < 0. Again, 
this is related to the finite distance and width of the injectors 
and can be reproduced in simulations. A similar asymmetry 
has been calculated in Ref. 17 (Fig. 12, the last s-wave plot) 
where the "inserted" flux is due to a trapped vortex, coupling 
<I>o/2 to the center of a conventional junction. If the size of 
the trapped vortex (in our case Aw and Ax) approaches zero, 
I C (H) approaches a perfectly symmetric shape with equal os- 
cillation periods at H > and H < (see Fig. 12 of Ref. 
the first d-wave plot). 

When a LJJ has a ^-discontinuity in its center, the ground 
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FIG. 5: Experimentally measured dependence +I C (H) for and 
+7r|0 and — 7r|0 states of the junction. The dependence of the neg- 
ative critical current —I C (H) is not shown, but it is perfectly mirror 
symmetric with respect to the if -axis. 
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FIG. 6: Sketch of a semifiuxon in a 0-7T-LJJ and two images outside 
the junction. Dashed line shows how the magnetic field of semi- 
fiuxon would look like in an infinite LJJ. Solid line shows the mag- 
netic field of a semifiuxon in a LJJ of L/\ j = 2. The arrow indicates 
the localized field produced by the injector. 



state consists of a semifiuxon spontaneously formed and 
pinned at the ^-discontinuity 1 71 1 1 lL 1 1 3W . as shown in Fig. [6] 
When L ~ A,/ the semifiuxon solution001, shown in Fig. [6] 
by a dashed line, is not valid anymore as it does not satisfy 
the boundary conditionsg^+L^) = 0. The semifiuxon so- 
lution derived in Refs. |7J, |g is valid in an infinite LJJ only. 
To construct the solution for a LJJ of finite length, we intro- 
duce images (anti-semifluxons) situated outside the left and 
the right edge of the junction at a distance of L/2 from each 
edge|28J. The final magnetic field profile 4> x (x) is shown in 
Fig-Elby the thick solid line. 

If one applies a bias current to a chain of antiferromagneti- 
cally ordered semifluxons, the semifluxons can rearrangel 11]: 
the positive semifiuxon at x = hops to the 7r-discontinuity 
at x — L and the image hops from x = L to x = 0. This pro- 
cess is accompanied by a transfer of flux equal to $n across 
the right edge of the junction. Subsequently the same mutual 
hopping happens at the left edge of the junction and the flux 
equal to <f>o is transferred across the left edge too. This is re- 
peated periodically in time, and the net transferred flux results 
in a finite voltage across the junction. This process is very 
similar to the dynamics which takes place in a LJJ biased at 
the zero-field step (ZFS)El. 

Presuming that the flux in our 0-7T-LJJ is transferred 
with the same maximum velocity, but only <I>o rather 2<£>o 
per hopping/reflection, we should observe half-integer ZFS 
on the IVC of the 0-7T-LJJ, as was recently predicted by 
Stefanakis llSII . Note that to observe semi-integer ZFS the 0- 
7T-LJJ should be not very long (L < 5A,/) so that the semi- 
fluxons can still hop over this distance. 

In Fig. we show two IVCs taken at H = 0: one with a 
classical ZFS recorded at 7j n j = and another with a semi- 
integer ZFS registered at /j n j = +13.5 mA. The asymptotic 
voltage of the ZFS is Vzfs = 388 /iV, while for the semi- 
integer ZFS it is ^3/2 = 196 ^V. Thus, 2F ZFS / 2 ~ V ZF s 
with accuracy of 1 %. This is the experimental confirma- 
tion of the existence of semi-integer ZFS 1 1 8 ] and semifiuxon 
hopping 1 11]. We also observe resonant structures at the half- 
integer ZFS, which is similar to structures we see in simula- 
tions. This kind of resonances typically appear when magnetic 



4 




-j— i— i— i— i— |— i— i— i— i— |— i— i— i— i— |— i— i 

100 200 300 

voltage across LJJ (^V) 



FIG. 7: Experimentally measured integer and semi-integer zero field 
steps on the IVC of the artificial 0-7T-LJJ at 7 inj = +13.5 mA (full 
circles) and Jj n j = (open circles). 

flux moves in a LJJ with inhomogeneities| 19]. 

We also measured the behavior of the semi-integer ZFS as 
a function of 7 in j and found that the asymptotic voltage of the 
semi-integer ZFS does not depend on /; n j . When k ^ 7r, one 
has two types of fractional vortices in the system: a vortex n 
at x = and two images (antivortices) k — 2ir at x = ±L. 
When the vortex and antivortex hop and exchange positions, 
the total flux transferred is still $o (they exchange a virtual 
fluxon| 11]). Although the voltage of the semi-integer ZFS is 
not affected by ij n j, the amplitude of the step I max changes 
with Jinj until the step gets completely hidden by the critical 
current. The dependence of the maximum current of the semi- 
integer ZFS on ijnj is shown in Fig. |4] No high order ZFSs 
were observed since our LJJ is rather short (~ 2A j). 

In conclusion, we have proposed, implemented and tested 
experimentally long Josephson junctions based on Nb-AlO^- 
Nb technology with an arbitrary ^-discontinuity of the Joseph- 
son phase created by passing a current through a pair of injec- 
tors. Formally this system is equivalent to a 0-ir long Joseph- 
son junction made of unconventional superconductors and al- 
lows one to study the dynamics of fractional vortices. We con- 
firmed the hopping of semifluxons under the action of applied 
bias current which manifests itself as the appearance of semi- 
integer zero field steps in the current- voltage characteristics of 
the junctions. 

Among the many interesting physical questions concern- 
ing the physics of fractional vortices we just mention two ex- 
amples for which arbitrary quantization and/or low damping 
are essential: (1) How does the eigen-frequency of a frac- 
tional vortex depend on the flux carried by it? How does the 
coupling between neighboring fractional vortices affect their 
eigen-modes (eigen-mode splitting) and can one construct an 
artificial ID crystal with controllable energy bands? (2) Does 
a semifluxon behave as a spin-i particle? Can it make quan- 
tum transitions f^l or be in a superposition of both states or 
does flux conservation prevent this? Can one create superpo- 
sitions of states such as Ti^iT f° r two coupled semifluxons? 



Finally, we note that an annular system with e.g. three semi- 
fluxons may serve as an interesting system to study frustration 
effects. 
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